Quantum critical transport in the unitary Fermi gas 



(N 

O 

(N 

1—5 

m 



C/2 



I 

o 



(N 
> 
O 
00 

o 



X 



Tilman Enss 

Physik Department, Technische Universitdt Miinchen, D-85747 Garching, Germany 

The thermodynamic and transport properties of the unitary Fermi gas at finite temperature T 
are governed by a quantum critical point at T = and zero density. We compute the universal 
shear viscosity to entropy ratio 77/s in the high-temperature quantum critical regime T 3> and 
find that this strongly coupled quantum fluid comes close to perfect fluidity 77/s = /i/(47rfc_B). Using 
a controlled large- A'^ expansion we show that already at the flrst non-trivial order the equation of 
state and the Tan contact density C agree well with the most recent experimental measurements 
and theoretical Luttinger-Ward and Bold Diagrammatic Monte Carlo calculations. 

PACS numbers: 03.75.Ss, 05.30.Fk, 51.20.-|-d 



I. INTRODUCTION 



The unitary Fermi gas is a basic many-body prob- 
lem which describes strongly interacting fermions rang- 
ing from ultracold atoms near a Feshbach resonance [l|- 
Q to dilute neutron matter. The properties in the dilute 
limit are independent of the microscopic details of the in- 
teraction potential and share a common universal phase 
diagram. A quantum critical point (QCP) at zero tem- 
perature governs the critical behavior in the whole phase 
diagram as a function of temperature T, chemical po- 
tential /i, detuning from the Feshbach resonance i/, and 
magnetic field h [4-6]. Whereas conventional QCPs sep- 
arate two phases of finite density, in our case the density 
itself is the order parameter which vanishes for ^ < 
and assumes a finite value for > Q. In the spin 
balanced case /i = 0, and at resonance v = Q the Fermi 
gas is unitary and scale invariant. In terms of the ther- 
mal length At = ft(27r/mfcsT)^/^ the density equation of 
state nXj, = /n(M/fcsT') is a universal function which has 
been measured experimentally [7.i lSj] ■ The unitary Fermi 
gas becomes superfluid at a universal Tc{^) « 0.4 /x Q, 
see Fig. [T] In this work we focus on the quantum criti- 
cal regime T > above the QCP &t h — v — {) and 



/Lt = 0, where nX 



3 _ 



/n(0) ~ 2.9 is a universal constant. 



Since the thermal length At is comparable to the mean 
particle spacing n~^/'^, quantum and thermal effects are 
equally important. There is no small parameter, and it 
is a theoretical challenge to compute the critical proper- 
ties. Recent measurements [8] and computations [13] 
of the equation of state now agree to the percent level. 
However, a precise determination of transport properties 
is much more demanding. 

In order to reliably estimate transport coefficients we 
perform controlled calculations in a large-A^ expansion 
@, Ell • Due to the lack of an intrinsic small parameter 
we introduce an artificial small parameter, which 
organizes the different diagrammatic contributions, or 
scattering processes, into orders of The original 

theory is recovered in the limit = 1. One can perform 
controlled calculations by including all diagrams up to a 
certain order in and these approximations can be 

systematically improved by going to higher order. This 
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FIG. 1: Universal phase diagram of the unitary Fermi gas. 



approach is similar to the e expansion in the dimension 
of space. The advantage over perturbation theory is that 
it is controlled even at strong interaction, while in con- 
trast to Quantum Monte Carlo it works directly in the 
thermodynamic limit and needs no finite size scaling. 

We thus obtain new results for the Tan contact den- 
sity [T^Jli'l and the transport properties in the quantum 
critical region. The shear viscosity ry = hX^'^ fridi/kBT) 
assumes a universal value at /i = 0. In kinetic theory 
r] = Pt is given by the pressure P times the viscous 
scattering time r, which is related to the incoherent re- 
laxation time of the gapless critical excitations above the 



QCP. The entropy density s = ks^T fs{^J■/kBT) at ^ = 



is exactly proportional to the pressure, s 
the viscosity to entropy ratio (at = 1) 



-Tt w 0.74 



kB 



5P/2T, and 



(1) 



is a universal number independent of temperature. A tem- 
perature independent ratio rj/s = h/^AnkB) has been 
found in certain string theories [1^1 and is conjectured 
to hold as a lower bound in other models Strongly 
interacting quantum fluids which saturate this bound are 
called perfect fluids [13] • Among real non-relativistic flu- 
ids the unitary Fermi gas comes closest to the bound and 
is almost perfect [T8-.20] , while for graphene the viscosity 
decreases logarithmically with temperature in the quan- 
tum critical regime [2l| . 
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Experiment Large-A'^ LuttWard BoldDiagMC 



2.674 3.108 [26] 2.90(5) [9] 

0.928 0.863 [26] 0.90(2) [9] 

2.320 2.177 [26] 2.25(5) [9] 

0.0789 0.084 [18] 0.080(5) [27] 



nAf, 2.966(35) [8] 

P [nksT] 0.891(19) [8] 
s [nfcs] 2.227(38) [8] 
C[k%] 

r]/s [h/kg] 1.0(2) [19, 28] 0.741 0.708 [18] 

TABLE I: Thermodynamic properties and transport coeffi- 
cients of the unitary Fermi gas in the quantum critical region 
= 0, r > 0: density n, pressure P, entropy density s. Tan 
contact density C, and shear viscosity rj, with Fermi momen- 
tum kp = {3n^n)^^^. Large- A'' results extrapolated to A'^ = 1. 



We compare our large- results at TV = 1 [22] with 
experimental measurements [l^ [2^ [13] and other the- 
oretical approaches, including self-consistent Luttinger- 
Ward [11 [p, nil and Bold Diagrammatic Monte Carlo 
(BDMC) [^calculations, see Table HI 

The excellent agreement between experiment and 
BDMC provides a reliable reference to assess the accu- 
racy of other methods. We find very good agreement 
of the pressure P with large- iV (3% above BDMC) and 
Luttinger-Ward (4% below) calculations, just slightly 
outside the error bars, and we find similarly good agree- 
ment for the entropy densitys. From the BDMC equa- 
tion of state simulations of (9|, one can extract (via the 
pair pro pagator) a preliminary value for the contact den- 
sity [23 C/kjr = 0.080(5). Our large-iV value is just 
1.4% below the BDMC value, which is remarkable given 
how simple the calculation is, while the Luttinger-Ward 
value lies about 5% above the BDMC value, just inside 
the error bars. Experimental measurements of the con- 
tact 24] yield C = 0.030(6) kp for the trapped gas at 
11 = {T/Tp = 0.64), which agrees well with trap aver- 
aged calculations [24|. However, knowledge of the trap 
averaged contact does not allow us to reconstruct the 
corresponding value for the homogeneous system, so we 
refrain from a direct comparison. Dynamical and trans- 
port properties such as 77/s are harder to compute than 
thermodynamic properties, which makes simple approxi- 
mations all the more valuable: we find that rj/s agrees to 
5% between large- iV and Luttinger-Ward theory, giving a 
narrow estimate. The viscosity of a trapped gas has been 
measured experimentally and agrees with trap averaged 
calculations but differs fr om the viscosity of 

the homogeneous system. 

The body of this paper explains how these values are 
obtained: in section we review the renormalization 
group (RG) analysis of the unitary Fermi gas and its 
universal phase diagram, in section IIIII we perform ther- 
modynamic and transport calculations using the con- 
trolled large- iV expansion, and in section IIVI we extract 
the /X = data from the self-consistent Luttinger-Ward 
calculation, before concluding in section|Vl In particular, 
in Appendix |^ we give a new derivation of the Tan adi- 
abatic and energy relations and show that they are sat- 
isfied exactly in self-consistent Luttinger-Ward approxi- 



mations, while Appendix IB] provides technical details on 
the quantum kinetic equation. 



II. PHASE DIAGRAM OF THE UNITARY 
FERMI GAS 

The interacting two-component Fermi gas is described 
by the action 



Sf = J d^x dr (a, _ Z_ _ 

goV'IV'I^i^t} (2) 



+ 



where ipa are Grassmann variables representing fermion 
species cr = i of equal mass m, and the imaginary time 
T = . . . /3 runs up to the inverse temperature /? = 1/T 
(we use units where h = 1 = ks)- Mo- is the chemical 
potential of species cr, but we will only consider the spin- 
balanced case /i = /J,^ = ^4,. 

In d = 3 dimensions the scattering amplitude for small 
relative momenta k can be written in the form [^ 



fik) 



1 



-1/a -ik + 



(3) 



where the scattering length a can be varied experimen- 
tally by an applied magnetic field, and the effective range 
Tg depends on the details of the interatomic potential. By 
fine-tuning to a Feshbach resonance 1/a — >■ the two- 
particle scattering remains strong at low energy A: — > 
and reaches the unitarity limit f{k) = i/k independent 
of Tg. The low-energy properties remain universal at fi- 
nite density rt > if re is much shorter than the mean 
particle spacing n~^/^. This condition kpre — >■ is real- 
ized physically for a dilute gas and near a broad Feshbach 
resonance as in ^Li [2]. 

A finite regularizes the contact interaction at short 
distances (UV), and for a sharp momentum cutoff A ^ 
l/]re] the detuning is related to the bare coupling go 
in ^ by 
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(4) 



Note that the resonance 1^ = can only be reached for 
attractive interactions ga < 0, when a bound state of the 
interatomic potential is at the continuum threshold. 

More generally, this can be understood from an RG 
analysis of the model ([2]): at zero temperature and den- 
sity the running coupling g obeys the exact flow equation 



m 



- = i2-d)g-- 



(5) 



which in 2 < d < 4 has an unstable fixed point at 
g^ = —2{d — 2) < corresponding to the Feshbach reso- 
nance. For smaller g < the fermions will form a BEG 
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of fermion pairs; for larger g > the flow runs toward 
the attractive fixed point 17 = of the free Fermi gas 
(BCS limit). At the Feshbach resonance fixed point the 
detuning v is a relevant perturbation with scaling dimen- 
sion dim[i'] = d — 2. 

The zero temperature phase diagram exhibits a quan- 
tum critical point at the Feshbach resonance ^ = 0, zero 
chemical potential /Lt = 0, and zero spin imbalance h = 0, 
where the "magnetic field" h couples to the difference 
in chemical potential — This critical point deter- 
mines a universal phase diagram for finite T, v, fi and h 
0, 0] . In this work we concentrate on the spin balanced 
gas h = at unitarity ^ = 0: the phase diagram for finite 
T and fi is depicted in Fig. [TJ 

On the lower right for n/T > {^/T)c there is a su- 
perfiuid phase of fermion pairs, while the left part is a 
normal Fermi liquid phase at finite density. The phase 
transition line Tc{^) ~ 0.4^ Q is universal and strictly 
linear, in contrast to the corresponding phase diagram 
for a dilute Bose gas On the left for /i/T — > —00 the 
Fermi liquid crosses over to a dilute classical gas. The 
line T = 0, /z < has zero density (vacuum). Here we 
focus on the high-temperature quantum critical regime 
T ^ and in the following we compute the thermody- 
namic and transport properties specifically for the repre- 
sentative value ^ = 0. 

It is useful to perform a Hubbard-Stratonovich trans- 
formation to decouple the fermion interaction. We intro- 
duce a complex field r) representing a fermion pair 
and write the Bose- Fermi action 

SbF = j d'^xdri^ r„ (^r-^~ f^a)^a 

(7 

_1|0|2_^^|^|_^*^^^^|. (6) 

Note that the pairing field has a positive gap because 
go < near the Feshbach resonance. The action ([5]) 
has the same critical behavior as the two-channel atom- 
molecule model at its zero-range fixed point 2, 5]. 

One can now proceed by integrating out the fermions 
to obtain an effective bosonic action for the pairing field 
(j). This action has bosonic vertices with any even number 
2n of fields which are given by a bare fermion loop with 
2n vertex insertions. In contrast to the repulsive Fermi 
gas, where these vertices are irrelevant in the RG sense, 
for the unitary Fermi gas these vertices all have marginal 
scaling. Already the particle-particle loop (n = 1), which 
contributes to the self-energy of the 4> field, changes the 
bare scaling dimension dim[0] = d/2 of the (p field by 
an anomalous contribution rj^ — A — d to the true scaling 
dimension dim[0] — {d + ri^)/2 — 2, which is independent 
of d (for 2 < d < 4). Similarly, all higher bosonic ver- 
tices n > 1 are singular for small external frequencies and 
momenta and scale marginally in the RG sense. There 
is no small parameter to suppress these higher loop di- 
agrams, and they are a priori equally important in the 
infrared (IR) . At zero density the 2n > 2 particle (3 func- 
tions are decoupled from the 2n = 2 particle /3 function 



in Eq. ([5]), which is therefore exact. Nevertheless, there 
may also be a three-particle resonance (Efimov effect) in 
the three-particle /3 function depending on the mass ra- 
tio and whether the particles are fermions or bosons [2^ . 
This changes the ground state from a two-particle to a 
three-particle bound state and leads to limit cycles in the 

RG flow fa^. 

At finite density all higher bosonic vertices couple back 
into the self-energy of the (j) field. In order to assess the 
quantitative importance of these higher vertices, one can 
introduce an artificial expansion parameter such as the 
dimension e = 4 — (ifor2<d< 4 Isif or 1/N for a 
large number of fermion flavors N Alternatively, 
one can use a Monte Carlo sampling of diagrams @ . In 
this work we perform a large- iV expansion and compare 
it with the results from other approaches. 

III. LARGE-N EXPANSION 

We modify the Bose- Fermi action 1^ by introducing N 
identical copies, or flavors, of f and ^ fermions, denoted 
by V'tra with a — and the flavor index a = 1, . . . , A^. 
The pairing field (j) is chosen to create an fi pair of any 
fiavor, and we obtain the action |Tl| 

SbF ^ I d'^xdri^^tpl^^dr - - lJ.a)lpaa 

era 

a a 

This action is 0{N) invariant under rotations in flavor 
space. The Gaussian integral over the fermion fleld yields 
the effective bosonic action 




.90 J 

= A^T^^{;^lnGo.(fc,w™) 

- T-\k,Um)mk,^n.)\^ + Om^^)} (8) 

with the trace running over the spin index a. The bare 
Fermi propagator Goo-(fc, aj,„) is given by 

Gg^ (fc, Wm) = ~iUJ„i + Ek - Ha (9) 

with dispersion Sk — k"^ /{2m), and the bosonic propa- 
gator —T{k,LUm) is given by the regularized T-matrix in 
medium, 

r"\fc,a;„0 = — +T V [ -^^Got(p,e„) 
go (27r)'' 

X Goi(k - PjWm - en) • (10) 
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The number of flavors N appears only as a global prefac- 
tor in the action ([5]), hence a controlled loop expansion 
is possible 5|. Each closed fermion loop contributes a 
factor of N, while each (j) propagator is suppressed by 
Even though the higher bosonic vertices still have 
marginal scaling, their contributions to the grand poten- 
tial are now suppressed quantitatively by powers of 
For T 3> /i the system is in the normal phase, and the ac- 
tion (HI) has a saddle point at {(j)) = 0. To order 0{1/N) 
the grand potential reads 



N 



r^^{21nGo(fc,c^™)-^lnr(fc,w„0} (H) 



Note that this order of the 1 /N expansion extrapolated 
to iV = 1 is exactly the Nozieres-Schmitt-Rink (NSR) 
theory [2^ . The Matsubara frequency summation can 
be continued analytically to real frequency. 



^|-2rin[l + e"''(^'=-^)] 

k 

/^^&H'5(fc,^,A^,i^)} (12) 

with the scattering phase shift 5{k,uj,^^v) — 
ImlnT(fc, w, /i, i^) and the Bose function h{uj) = 
[exp(/?a;) — Specifically in d = 3 the T-matrix reads 
(in the spin-balanced case h = Q) 



T-\Kuj) = - 



mv 
An 
m 
27r2fc 

X In 



^3/2 

47r 



— cj — iO — 2/i 



P 



1 -(_ g/3(ep-M) 
iO + 2/X - £p - Sk-p 

Lj + iO + 2/1- Ep - Sk+p 





LJ - 



(13) 



The integral is convergent and readily evaluated numer- 
ically. 



A. Thermodynamics 

Using we obtain for the pressure P ~ —Q./L'^ 
(equation of state) at ^ ^ Q, v — h — Q and T > 0: 

Z = ^ = p(o) + 1 p(i) + . . . 

N NL<^ N 
1 



where 



= (1.734 400+ —0.747 561 jTA^;'' (14) 



p(o) = 2(1 - 2-='/2)C(5/2) TA^3 



Since the unitary Fermi gas is scale invariant the internal 
energy density e is proportional to the pressure (3^ 

£ 3P / 1 \ 

— = = (2.601600+ —1.121341 rA;;^ (15) 

Also the entropy density s ^ dP/dT = (e + P - ^n) /T 
at unitarity and /i = is proportional to the pressure, 

s 5P / 1 \ 

— = = (4.335 999+— 1.868 902) A-3. (i6) 

The density at ^ = to order 0{l/N) is 

1 



^ ^ d{P/N) ^ 
N dji 

1 

N' 



,(1) 



N 



1.530 294+ ^1.143 936) A;^^ 



(17) 



where 



=2(1- 2-1/2)^(3/2) 
/IN f d^k duj , , .d6(k,uj) 

j j2nf-'^''^^r- 

If this order of the 1 /N expansion is evaluated at iV = 1 
(NSR) we obtain for the density 



71 = 2.674 230 A^;^ (^^ = 1) 



(18) 



The ratio of thermal length to mean particle spacing, 
Xtti^^^ ~ 1.388, is of order unity, hence quantum and 
thermal fluctuations are equally important in the high- 
temperature quantum critical region. The density deter- 
mines the Fermi temperature 



ksTp = — 
2m 



(37r2n) 



2ra 



(19) 



which is useful to compare with data given in terms of 
the reduced temperature 

T /Sa/tt \ -2/3 
e= — ^[^n\l) = 0.681496 (TV = 1). (20) 

Finally, the Tan contact density is defined as the total 
spectral weight (density) of the pairing field (T2| - [l^ 



C — m 



2/J,*, 



m2 f d^k duj 

— — &(a;) lmT(K, uj) 



N J (27r)3 TT 



26.840128^ 

N 



(21) 



At = 1 the contact can be expressed in terms of kp 
using Eq. (HH) which yields C = 0.0789 k%. This is equiv- 
alent to the non-self-consistent T matrix result [s^l and 
agrees with the BDMC calculation within 1.4% (see Ta- 
ble |T]), but it differs from the result in [6i] by a factor of 
two. 

Note that the Tan adiabatic theorem llSl 



d{-P/N) C 



dv 



47rm 



(22) 
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is fulfilled exactly in the \/N expansion: the change of 
the pressure with detuning is 

d{-P/N) _ I r (fk duj dS{k,u}) 

^ J (27r)3 V 

m f d^k duj ^ , . , 

b{uj)lmTik,u}) (23) 



Attn J (27r)3 tt 



because the change of scattering phase shift with de- 
tuning is d5{k,uj) / dv = (m/47r) IniT(fc, w), and using 
Eq. dUl) we obtain (j^ . 



B. Transport 

At iV = cx) the fermions are free: once a shear flow is 
excited in the infinite system it will continue forever, and 
the dynamic shear viscosity is 



ri{uj) = TrPS{uj) . 



(24) 



The Drude weight is proportional to the pressure, in ac- 
cordance with the viscosity sum rule [18., 34, 35|. At order 
the fermions acquire a self-energy correction by scat- 
tering off pairing fluctuations, so for large N the fermions 
are almost free quasi-particles with lifetime 0{N) and an 
energy shift of the quasi-particle dispersion ReS ~ 
In kinetic theory the dynamic viscosity becomes 



Pt 



1 + [uJTf 



(25) 



with the viscous scattering time r = 0{N): the 5{uj) 
function in p4|) is broadened to a peak of width 
and height N . Note that the high-frequencytail r] ^ 
C/lbn^ymu; is not seen in kinetic theorv [3^ . 

In order to compute transport properties for large N 
it is justified to use the quantum Boltzmann equation 
0, mi '■ (i) the fermions propagate as free particles be- 
tween collisions, up to subleading corrections, (ii) the 
collision integral ImS ~ 1/iV contains only particle- 
particle scattering described by the medium T-matrix 
T{k,uj) [Eq. (fTO|) ] because particle-hole scattering ap- 
pears at higher orders, and (iii) in addition to the col- 
lision (dynamic) term there is a shift of the dispersion 
(kinetic) term ReS ^ 1/A^ of the same order. However, 
it is only a subleading correction to the leading real term 
^xy/N ^ (see below) and can be neglected. 

Based on these considerations we arrive at the Boltz- 
mann equation [3^ [s^l 

df ^. df df 1 

for the distribution function /(p,r,i), where /[/] is the 
collision integral. For the shear viscosity we consider a 
velocity field u = Ux{y)^ with a small shear gradient 
dux/dy, and the local equilibrium distribution /(p) = 



/"(e — u • p) with e = /2m. In the stationary limit the 
Boltzmann equation becomes [13] 



(27) 



df 1 



The velocity gradient induces a momentum current den- 
sity 

f d^n 8ii 
n.. = 2iVy (2^-.P./(p) = -'7^ (28) 

proportional to du^/dy, with the coefficient given by the 
shear viscosity rj. We choose a deviation from the equilib- 
rium distribution, f = f + Sf with Sf = f{l ~ f)(p{p) 
and ip{p) — VyPx/T, such that the momentum current 
density is 



n. 



2N 



(27r)3 v^^-'P ^ ■'P' 



= P. 



(29) 



This is equal to the pressure for free fermions (N = oo) at 
arbitrary temperature, as can be seen by integrating by 
parts. We can now replace —dux/dy — P/r] in (P?]) and 
take moments of the Boltzmann equation by integrating 
both sides with 2N J d^p/ {2itY VyPx- The left-hand side 
becomes 



2NP 



d^p 
TjT J (27r)3 



p2 



(30) 



while the right-hand side yields the collision integral |35l . 

113,131 



Cxy 2 



d^p 



(27r)3 

2 r d^p 

T 



VyPxI[Sf] 



(2^)3 ''^^^ 7 (27r)3 
x/X(l^/p')(l-/°i) 

X [lp{p) + Lp{pi) ~ Lp{p') - if{p[)] 



d^pi 



(31) 



where fermions with incoming momenta p, pi scatter 
into outgoing momenta p', p[. It will be convenient to 
express these momenta in terms of the total momentum 
q — p + Pi and the relative momenta k = (p — Pi)/2 
(k' = (p' — pi)/2) of the incoming (outgoing) particles, 
with |k'| = |k| by energy conservation. The occupation 
numbers give the probability that the incoming states are 
occupied, and the outgoing states are not. The differen- 
tial cross section is given by the medium T-matrix 



da 

dh 



in 



Tip + Pi, UJ = Sp + Ep^ - 2y) 



(32) 



In the vacuum limit the center-of-mass scattering de- 
pends only on the relative momentum fc. 



da or 
dh ^ TToFk^ 



(vacuum) 



(33) 
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but at finite density there is an additional dependence on 
the total momentum q in the medium T-matrix T{q, k) — 
T{q, uj — 2eq/2 + 2sk — 2m)- In relative coordinates the 
shear term in Eq. ([5T|) is [39| 



The collision integral then reads 
2 f (Pq f cPk 



k'^k'y 



mT/2 



(34) 



" nniT 



(27r)3 J (27r)3^'^^^''''^''^''^'''^-^'i/2+k/q/2-k 



dflk 



An 



-{kxky - k'^ky){l /q/2+k')(l ~ /q/2-k') 



1 



dqq^ / dkk'^\T{q,k)\'' 



(35) 



with the i?-wave angular average Ie{q,k) over the Fermi 
distribution functions derived analytically in appendix IbI 
(the d-wave average ^^1^2(9'^) contributes only 0.2% to 



the integral (|35|) ). Finally, only two integrals over the 
radial momenta q and k have to be performed. In the 
dilute classical regime the collision integral can be com- 
puted analytically, 



32^/2z2^2A^^ 



(36) 



with fugacity z — exp(/3/i), and in the same limit the 
pressure is Pel = 2zTX^^N. The viscosity is then given 
by M 



^ xy 



8%/2 



4.165 203 A^^A^^ 



(37) 



In the high-temperature quantum critical regime T > 0, 
/i = the collision integral has to be computed with the 
full medium T-matrix T{q,uj) from Eq. (|10p . which is 
done numerically and yields 



0.935 683 T^Ay^ 



(38) 



and together with the pressure at leading order in 1/iV, 
P = 7VP(°) = 1.734 400 TA^^A^, we obtain in the quan- 
tum critical regime 



Cxy 



3.214 917 Ay^A^^ 



(39) 



This value is about 20% lower than in the dilute classical 
limit ([57)) . which is mostly due to the reduced pressure, 
while the effects of reduced density and increased medium 
scattering almost cancel each other in Cxy ■ With the vis- 
cous relaxation time t = P/Cxy and the entropy density 
s = 5P/(2r) we obtain the universal viscosity to entropy 
ratio independent of temperature. 



0.741448- 



hN 



(40) 



A related computation of the viscosity using the medium 
T-matrix has been performed for large attractive inter- 
action kpa = —11.8 which found f] = 2.3Hn for ^ = 
at T/Tf = 0.7 fsl, slightly larger than our value ^ at 
= 1 . Note that we have evaluated rj using only a single 
moment of the Boltzmann equation ()27p , but is has been 
shown that corrections to rj from higher moments are less 
than 2% [i^. A similar transport calculation using the 
medium T-matrix in two dimensions has been performed 
recently (4ll |. 



IV. LUTTINGER-WARD THEORY 

The Luttinger-Ward theory provides a systematic way 
to obtain self-consistent and conserving approximations, 
such that the Green's functions satisfy all symmetries and 
conservation laws of the model [i^, ll^ . The Luttinger- 
Ward functional ^[Ga, Gb] can be defined in terms of full 
fermionic propagators Go- and full bosonic propagators 
Gb = —T- The exact theory is given by an infinite 
set of irreducible contributions to the $ functional which 
cannot be evaluated in practice, so typically one chooses 
a subclass of diagrams. For the unitary Fermi gas a very 
successful approximation is to use ladder diagrams with 
full fermionic Green's functions [1^, [2^ . Then the full 
T-matrix is given by an expression similar to (jlOp but 
with full Green's functions, 

r-\k,.^n)-y^+Tj:J^G,{p,e„) 

X G^(k- p,w„ - e„) . (41) 

Since we are interested in the high-temperature critical 
region we consider only the expressions valid in the nor- 
mal phase. The Luttinger-Ward theory then prescribes 
that the f fermionic self-energy is given by scattering a 
4, fermion off pair fluctuations described by the full T- 
matrix, 

(42) 

and analogously for E^. The Dyson equation determines 
the full fermionic Green's functions 



G^^{k,UJri 



(43) 



This set of equations (j4T |) -(|43 l) is solved self-consistently 
by iteration [25, 26]. The resulting Green's functions in 
Matsubara frequency can be continued analytically to ob- 
tain the spectral functions in real frequency, which show 
substantial broadening near T^, and additional excitations 
beyond a single quasi-particle peak f4^-|46|. Similar fea- 
tures are observed in the spin polarized case [47l - l50| . 
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The pressure P 
potential [26j 



-Vl/L'^ is obtained from the grand 



r! = T^^{^hiG,(fc,c<.„0 



^[1 - G^^{k,Lo,^)G„{k,Lo^)] - lnr(A;,c^„,)} (44) 



Acknowledgments 

I wish to thank Lars Fritz, Subir Sachdev, Jorg 
Schmahan, Richard Schmidt, and Wilhelm Zwerger for 
fruitful discussions and Mark Ku, Thomas Schafer, Chris 
Vale, Felix Werner, and Martin Zwierlein for sharing 
their data. 



evaluated using the self-consistent fermion propagator 
and the full T-matrix. We extract the high-temperature 
quantum critical behavior from the existing thermody- 
namic data pGl j interpolated at = 0. Specifically, we 
make a cubic spline interpolation of /i(T) and find the so- 
lution of ^J,{T) = at 6* = T/Tp = 0.6165, which implies 
nX^ = 8/(30r)6'-3/2 = 3JO8. Furthermore, we find 
P = 0.8630 n/cijT, s = 2.177 nks, and C = 0.084 353 fc|,, 
which can be recast in terms of At. These values are 
summarized in Table U and are remarkably close to the 
experimental values. 

The shear viscosity ri{T,uj) has been computed in 
Luttinger-Ward theory as a function of temperature and 
frequency 18| : it has a Lorentzian peak at low frequency, 
followed by a universal tail ri{T,uj) ~ C{T)/15w^/rriuj 
proportional to the contact density. We make a cu- 
bic spline interpolation of = and find the root 
at ri{T,uj = 0) = 1.5409 which yields 77/s = 
0.7077 h/ks- This result is slightly lower than the large- 
TV value in Eq. (|40p. We note that in this self-consistent 
calculation the minimum of 77/s O.&h/kB is found at a 
somewhat lower temperature T/Tp ss 0.4 fisj . 



V. CONCLUSIONS 

The unitary Fermi gas in the high-temperature quan- 
tum critical region is a challenging many-body problem. 
It is strongly interacting, with the density almost twice 
the non-interacting value at /x = Q, and has no small 
expansion parameter. Still, our large-A^ results at the 
first non-trivial order beyond the free Fermi gas are al- 
ready remarkably close to reliable experimental and theo- 
retical results 8, 3\ ■ A main result of the present paper is 
that this is true also for the transport properties rj/s once 
medium effects are included in the quantum kinetic equa- 
tion. A possible reason for this good agreement is that 
large- A'^ and Luttinger-Ward approximations satisfy the 
Tan adiabatic and energy relations exactly, as we show in 
Appendix [X] In addition, Luttinger-Ward theory exactly 
fulfills the scale invariance of the unitary Fermi gas (181 ). 
For a better comparison between calculations for the ho- 
mogeneous system and experiments it would be desirable 
to have local measurements in the spirit of Ref. Q also 
for the contact and transport properties, since the com- 
parison of trap averaged quantities is less sensitive to the 
details of the temperature dependence. A promising step 
in this direction is to selectively probe atoms near the 
center of the trap in order to extract the contact density 
from the tail of the momentum distribution 1511 . 



Appendix A: Exact Tan relations in Luttinger-Ward 
theory 

Consider the fermionic action Q : a small variation of 



the quadratic term, 5Gq 
grand potential 



will lead to a change in the 



= -tr(G(5Go^) 



(Al) 



with the trace running over space, time and possibly spin 
indices. However, this equation is often violated if ap- 
proximations are made for the full Green's function G. 
A unique feature of conserving approximations, which are 
derived from a Luttinger-Ward functional ^[G], is that 
Eq^ (|A1[) holds exactly even for approximate Q, and G 

For the strongly interacting Fermi gas it is conve- 
nient to start from the Bose-Fermi action ^ and de- 
fine a Luttinger-Ward functional ^[GcGb] in terms of 
both fermionic and bosonic Green's functions jlSj, |25„ ^26)] . 
Then a variation of the microscopic parameters SG^^ 
and/or 6G^g induces a change of the grand potential 

M 



SQ^- tviGaSGo^) + tr(GB<5Goi) 



(A2) 



Again, this exact equation continues to hold within con- 
serving approximations with full self-consistent propaga- 
tors Go- and Gb- 

We will now show that the Tan adiabatic theorem llSll 



C 



d{-l/a) 47rm 
and the Tan energy formula [T2j 

^ f d^k / C 



k^ 



C 



4TTma 



(A3) 



(A4) 



are consequences of (|A2|1 and therefore hold not only in 
the exact theory but in any conserving approximation, 
including the self-consistent T-matrix approximation in- 
troduced in section ITVl A variation of detuning changes 
only the bosonic quadratic term 



^-iM N 1 m / 1 2A\ 
Gob A:,^,n = = ;p — + — 



(A5) 



in the action ([5]), 

9GoJ(fc,tj„) 
di-l/a) 



= 



dGg^{k,uj„,) _ m 
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The variation of the grand potential is then 

'd{-l/a) 



d{-l/a) 



= -tr 



Gf 



d{-l/a) 



in 



tr(GB) 



(A7) 



with the density of bosons expressed by the Tan contact 
density, 



Gb{x = Q,t^-Q) = {c^*4>) 



c 



(A8) 



Inserting (|A8|) into (|A7|) directly yields the adiabatic the- 
orem (|A3|) . In order to derive the energy formula we 
consider a variation of mass, 



.1 dn/L^ 
d{m^^ 



(A9) 



Usually this yields only the kinetic energy (cf. Eq. (61) 
in [Hlj), but in our case the interaction term Aira/m also 
depends on mass, so (jA9l) is the full internal energy e — 
(H) including the potential term. Specifically, 



m 



Sk 



-1 9GQ^{k,uj„ 

^i dG^g{k,uJm) _ _ 2A \ 

dni^^ 47r Va tt / 



(AlO) 
(All) 



and with the momentum distribution function 
— G„{k,uj„i) — rii^a- we obtain the internal 
energy density 



m 



L3 



E 



d^k 



m /I 2A\ G 



47r (( 



(A12) 



where the k integral extends to the momentum cutoff A. 
The regularization term AC/(27r^m) can be written as 
Scr / d^ k / {2tt)^ EkC / k'^ , and we arrive at the energy for- 
mula (jA4[) . In a similar way the Tan pressure relation 
has been derived in the Luttinger-Ward theory by an in- 
finitesimal scale transformation on the grand potential 
[m . This concludes our proof that the Tan relations are 
fulfilled exactly in the self-consistent T-matrix approxi- 
mation. 



Appendix B: Quantum kinetic equation 

A useful feature of d = 3 dimensions is that the angular 
averages of the distribution functions can be performed 
analytically. One can write the product of Fermi func- 
tions in (l35l) with Ik' I = Ikl as 



/(eq/2+k)/(eq/2-k)[l " / (eq/2+k' )] [1 " /(£q/2-k')] 

_ }_ (Til] 

4(cosha + cosh ^.t) (cosh a + cosh bx') 

with X = k • q, .t' = k' • q and a = (£q/2 + £k ^ m)/^: 
b = kq/ {2mT). The angular average over the solid angles 
of the vectors q, k and k' is then 

= '^{lLo{q,k)~ll,iq,k)) (B2) 



where we have defined the i'-wave angular average of the 
distribution functions 



1 I 

h{q,k) = - I dx — - 

4 J_i cosh a 



cosh bx 



(B3) 



with Legendre polynomials Pi{x). The s-wave average is 
given by 

T , i\ 1 1 cosh[(a + &)/2] 

'^-°('^'')=25^^ %osh[(a^5)/2] 

while the c?-wave average can be expressed in terms of 
polylogarithms Lis(z), 



Ii=2iq,k) = Ii=o{q;k) 



1 



46-^ sinh a 

6Li3(-e'^+'') - 6Li3(-e^-'^) - 66 Li2(-e'^+^) 

+ 66Li2(-e''"") -I- a{a^ ~ Sab + tt^)] . (B5) 



Thus, all angular integrations can be done analytically 
and only the two radial integrations over q and k in 
Eq. (|35p need to be performed numerically. 
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